Skyrmion Lattice in a Doped Semiconductor 
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We report a comprehensive small angle neutron scattering study (SANS) of the magnetic phase 
diagram of the doped semiconductor Fei-ccCoa^Si for x = 0.2 and 0.25. For magnetic field parallel 
to the neutron beam we observe a six-fold intensity pattern under field-cooling, which identifies the 
A-phase of Fei-ccCoxSi as a skyrmion lattice. The regime of the skyrmion lattice is highly hysteretic 
and extents over a wide temperature range, consistent with the site disorder of the Fe and Co atoms. 
Our study identifies Fei-ccCoxSi is a second material after MnSi in which a skyrmion lattice forms 
and establishes that skyrmion lattices may also occur in strongly doped semiconductors. 



PACS numbers: 72.80.Ga, 72.15.-v, 72.25.-b, 75.30.-m 

Recently a skyrmion lattice was identified in the cubic 
B20 system MnSi [Hll], that is, magnetic order represent- 
ing a crystallization of topologically stable, particle-like 
knots of the spin structure originally anticipated to occur 
in anisotropic materials [3]. This raises the question for 
further magnetic materials with skyrmion lattices and if 
they are a general phenomenon in cubic magnets without 
inversion symmetry as suggested by our theoretical treat- 
ment in [1 . Because MnSi is a pure metal, an additional 
question concerns if skyrmion lattices are sensitive to dis- 
order and whether they also exist in semiconductors and 
insulators. More generally, the microscopic identification 
of a skyrmion lattice in MnSi represents also a showcase 
for similar lattice structures considered in nuclear physics 
[H |5] , quantum Hall systems [6l [7] , and liquid crystals [8 . 

The B20 transition metal silicides TSi (T=Fe,Co,Mn) 
are ideally suited to clarify these questions. While FeSi 
is a nonmagnetic insulator with strong electronic corre- 
lations [9] its sibling CoSi is a diamagnetic metal [10]. 
For increasing Co-content the series Fei_a;Coa;Si displays 
a insulator to metal transition at x = 0.02 becoming in- 
creasingly metallic with a dome of helimagnetic order in 
the range 0.05 < x < 0.7 [11 . Being a strongly doped Si- 
based semiconductor the magneto-transport properties of 
Fei_a^Coa^Si attract great interest [T2j[T3]. 

The helimagnetic order in Fei-^^Coa^Si results from a 
hierarchy of energy scales, of which ferromagnetic ex- 
change on the strongest scale favors parallel spin align- 
ment and Dzyaloshinsky-Moriya spin-orbit interactions, 
permitted in the non-centrosymmetric B20 structure, fa- 
vor perpendicular spin alignment on a weaker scale. With 
increasing x the wavelength of the resulting helical mod- 
ulation varies from about 200 A to 2000 A dHHl]. The 
propagation direction of the helical modulation is deter- 
mined by crystal- field interactions on the weakest scale. 
Based on SANS experiments [Ml [15] |T6 and Lorentz 



force microscopy [17] it was concluded that the propa- 
gation axis of the helix in Fei-^^Coa^Si is (100) for all x. 
While this differs distinctly from the propagation direc- 
tion in MnSi being (111), it is perfectly consistent with 
the cubic crystal structure [I8l [19]. Further, while the 
DM interaction and the crystal structure have the same 
chirality in MnSi [20], they have opposite chirality in 
Fei_a^Coa;Si [21]. Finally, due to the strong doping in 
Fei_a;Coa;Si there must be also site disorder of the Fe 
and Co atoms even in excellent single crystals. Thus, 
although Fei_ajCoa;Si and MnSi exhibit rather similar 
magnetic properties, there are strong differences of the 
underlying electronic structure. 

The magnetic phase diagram of Fei-a^Co^^Si exhibits 
three prominent features [14l |15l [16]. First, a state with 
a helimagnetic modulation up to a critical field re- 
ferred to in the following as zero-field cooled (ZFC) state. 
Second, for a magnetic field in the range Bd < B < Bc2 
a helical modulation parallel to the magnetic field is sta- 
bilized, forming the so-called conical state. Third, in the 
so-called A-phase, a small pocket just below Tc, the heli- 
cal modulation is perpendicular to the applied magnetic 
field. Prior to the work reported here all studies were car- 
ried out with the field perpendicular to the neutron beam. 
Because this configuration is mostly sensitive to helical 
components parallel to the field and only selected helical 
components perpendicular to the field the information on 
the A-phase was incomplete and it was believed that the 
A-phase represents a single-Q helical state |T4 l [T5 l [T6]. 

Here we report a comprehensive SANS study of the 
magnetic phase diagram of Fei-a^Coa^Si for x = 0.2 
and X = 0.25. In contrast to previous SANS studies 
[m [151 [E] we explore the magnetic phase diagram also 
with the magnetic field parallel to the incident neutron 
beam, a configuration that is sensitive to all helical mod- 
ulations perpendicular to the field. We also determined 



2 




-0.02 0.00 0.02 . -0.02 0.00 0.02 

Qx(A") Q,(A-^) 



FIG. 1: Typical SANS data of Fei_^Co^Si for x = 0.2. (A) 
through (C) Intensity after zero field cooling (ZFC). (D) A- 
phase after ZFC for increasing B parallel to the neutron beam. 
(E) through (G) Six-fold intensity distribution in the A-phase 
after field-cooling for B parallel to the neutron beam. (H) A- 
phase after field-cooling for B perpendicular to the neutron 
beam. In all panels the intensity on the left-hand side is 
reduced due to a small inefficiency of the neutron detector. 

differences of the phase diagram for zero-field cooling 
(ZFC) and field-cooling (FC). As our main result we iden- 
tify the A-phase of Fei-a^Coa^Si as a hexagonal skyrmion 
lattice akin that seen in MnSi. By comparison to MnSi, 
the skyrmion lattice yields hysteretic features and may 
be observed over a large temperature range. Thus we 
identify a second example of a skyrmion lattice, how- 
ever, in a doped semiconductor with strong disorder. In 
addition we find that the ZFC state, long known to be 
metastable [15], is characterized by scattering intensity 
on the surface of a small sphere in reciprocal space with 
broad maxima in the (110) direction, reminiscent of the 
partial magnetic order in MnSi at high pressure [22] . 

For our study several single crystals were grown at 
TUM by optical float-zoning with a UHV-compatible im- 
age furnace [23 . Laue x-ray diffraction, EDX, polar- 
ized light microscopy and single-crystal neutron diffrac- 
tion established an excellent sample quality. The magne- 
totransport properties, magnetization and specific heat 
were in excellent agreement with the literature [13 as 
well as small single-crystals grown by vapor transport 
at MPI-CPfS, Dresden. All samples were carefully ori- 
ented using Laue x-ray diffraction. A large single crystal 



with X = 0.2, studied most extensively, was additionally 
oriented on the diffractometer RESI at FRM II. Using 
neutrons with a wavelength A = 1.0402 A ± 1% we deter- 
mined ten Bragg peaks that were fitted simultaneously. 
The sample had a lattice constant a = 4.483 ± 0.02 A 
with an excellent resolution-limited mosaic spread better 
than 0.2°. We further confirmed within the resolution 
limit that there were no superstructures, structural short- 
range order, or lattice distortions. The modulus of the 
helical modulation Q ^ 0.017 A~^, transition tempera- 
ture Tc ^ 30 K and the lattice constant were in excellent 
agreement with Refs.[T5l [T6[ |24], but differed from crys- 
tals grown in a tri-arc furnace p!4] . 

Our SANS studies were carried out at the diffractome- 
ter MIRA at FRM II, using neutrons with a wavelength 
A = 9.7A±5%. Backgrounds were determined at high 
temperatures and subtracted accordingly. A Cd marker 
on the sample support confirmed that the sample was 
always oriented correctly. Unfortunately the ^He delay- 
line detector displayed an inefficiency on the left-hand 
side, consistently visible in Fig.[l] This does not affect 
the conclusions of the work reported here. The sample 
was cooled with a pulse-tube cooler. Magnetic fields up 
to 0.3 T were applied with a bespoke set of Helmholtz- 
coils. Data were recorded for a fixed sample orientation 
following rocking scans with respect to a vertical axis of 
typically ±15°. 

In the following we focus on the ZFC state and the 
A-phase. Data recorded in the conical state are not dis- 
cussed in detail since they agree very well with previ- 
ous studies. To determine the three-dimensional distri- 
bution of scattering intensity of the ZFC state we have 
recorded data for different sample orientations. As shown 
in Fig.[l]A, B and C for (111), (110) and (100), the ZFC 
intensity for x = 0.2 is characterized by a uniform ring, 
a ring with minima for (100) (horizontal direction), and 
a ring with minima for (100) (diagonals), respectively. 
Taking together Figs.[l]A, B and C, the ZFC state is char- 
acterized by intensity on the surface of a small sphere in 
reciprocal space with broad intensity maxima for (110) 
[25^. This contrasts previous studies, where a similar 
ZFC intensity distribution was reported with maxima for 
(100) [m |15l |2^. We have therefore confirmed our 
results for several samples with x = 0.2 from different 
growths. Moreover, for x = 0.25 we find a similar scat- 
tering distribution in the ZFC state (data not shown). 
Thus data in Ref. [HI |16] may have either been indexed 
incorrectly or an unexplained sample dependence occurs 
for X = 0.2 and 0.25, where our samples by all accounts 
seem excellent. The broad intensity maxima for (110) in 
the ZFC state are also unusual, because the propagation 
direction of the helical modulation in the cubic B20 crys- 
tal structure in leading order may be either parallel to 
(111) or (100) [I8l[19]. In fact, the ZFC intensity distri- 
bution corresponds to the partially ordered state of MnSi 
[22] . which inspired proposals of spontaneously forming 
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FIG. 2: Magnetic phase diagram of Fei-a^Co^cSi for x = 0.2. 
Gray arrows show the location of temperature and field scans. 
Red data points show transition fields and temperatures as 
determined in these scans. Panel (A) through (C) show the 
magnetic phase diagrams after ZFC; in the green shaded area 
the scattering intensity of the ZFC state is observed. Panel 
(D) through (F) show the magnetic phase diagrams for FC. 
Across the dark blue area the intensity displays a maximum. 



spin structures such as skyrmion lattices [271 ESI ESl [30] . 

Shown in Fig.[l]D through G are typical intensity pat- 
terns in the A-phase. These patterns are exclusively seen 
in the plane perpendicular to the magnetic field. After 
ZFC and for increasing field we observe a uniform ring 
of intensity regardless of the crystallographic orientation 
(Fig.[l]D). In contrast, for FC the intensity pattern in the 
A-phase is characterized by a six-fold symmetry. For the 
(111) and (110) directions six spots are observed as shown 
in Fig.[l^ and F, respectively, where two of the spots are 
aligned with a (110) axis. Yet, for the (100) direction the 
six-fold pattern is composed of twelve spots comprising 
two sets of six spots, as shown in Fig.[Tp. This is char- 
acteristic of two domain populations, where each domain 
yields the same hexagonal symmetry seen for (111) and 
(110). Interestingly, two spots of each domain population 
are now parallel to (100). We have finally also confirmed 
the formation of the conical phase, which coexists with 
the A-phase in Fig.[l]H (horizontal and vertical spots are 
due to the A-phase and conical phase, respectively). 

Shown in Fig. [2] is the magnetic phase diagram for 
X = 0.2 (similar diagrams for x = 0.25 are not shown). 
We find excellent agreement with published work where 
data are available. Starting from a ZFC state we ob- 
serve: (i) scattering intensity with broad maxima for 
(110) up to Bel, where Bci{T 0) ^ 95 mT and 
dBci/dT ^ — 3mTK~^ are isotropic (green shading), 
(ii) a conical phase as reported before (gray shading), 
and (iii) a ring of intensity in the A-phase, where the 




.0 


B <111> ■ 


• 24 mT 




A31 mT 




▼ 37 mT 




044mT 




451 mT 




►58 mT 










5 10 


15 20 25 30 




T(K) 



20 40 60 80 
B(mT) 




10 15 20 25 30 
T(K) 



FIG. 3: Temperature and magnetic field dependences of the 
integrated scattering intensity. (A) Intensity for B \\ (110) 
and B perpendicular to the neutron beam. (B) Intensity for 
B II (111) and B parallel to the neutron beam. (C) and 
(D) Intensity for B parallel to neutron beam under FC. All 
panels use the same arbitrary scale; gray shading shows the 
background. 



wave vector is perpendicular to the field (blue shading). 
For FC the salient features are: (iv) the ZFC scattering 
intensity is never observed, (v) scattering intensity in the 
A-phase with a six-fold symmetry and wave vectors per- 
pendicular to the field (blue shading), and (vi) a small 
field range in which the A-phase exists down to the low- 
est temperatures (this is not seen for field in the (111) 
direction) . The regime where the intensity in the A-phase 
displays a maximum is shown as dark blue shading, while 
the remaining regime is shown in light blue shading. 

Metastable aspects of the ZFC state and the A-phase 
are best illustrated by the integrated intensity shown 
m Fig.[3| For field scans with B perpendicular to the 
neutron beam we reproduced published results [M]. As 
shown in Fig.|3]A for (110), there is essentially no inten- 
sity up to Bel for increasing B and low T above which 
the conical phase appears. In contrast, for decreasing 
field the intensity remains unchanged high below Bd . In 
the vicinity of Tc the intensity displays a reversible min- 
imum in the A-phase between Bai and Ba2^ where the 
intensity includes contributions of the A-phase and the 
conical phase (see Fig.[l]H). 

Scans for field parallel to the neutron beam shown in 
Fig.[3|B, C and D are the counterpart to that shown in 
Fig.[3|A, notably they show all helical components per- 
pendicular to the field. For (111) and increasing field 
the ZFC intensity well below Tc is unchanged until it 
vanishes above Bd (Fig.|3]B). However, in field scans 
of decreasing strength the ZFC intensity is not recov- 
ered (data not shown). Near Tc the intensity reversibly 
shows a maximum in the A-phase (data at 28 K). In tem- 
perature scans a corresponding maximum is seen in the 
A-phase (see Fig.[3]C for (111)). An additional feature 
emerges in temperature scans for (100) and (110), shown 
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in Fig.|3]D for (100). For 5 > 31mT and PC the inten- 
sity displays not only a maximum, but also an additional 
increase down to the lowest temperatures. The regime of 
the maximum near Tc agrees between ZFC and FC for all 
three directions (compare dark blue shading in Fig.|2]A, 
B and C with D, E and F, respectively). Interestingly, 
for (100) the A-phase even borders directly on the ZFC 
state (Fig.[2]C). 

The scattering pattern and hysteretic behavior re- 
ported here are fully consistent with a scenario where 
the A-phase is thermodynamically stable in the dark blue 
areas and metastable in the light blue region. The char- 
acterizing 6-fold pattern is, however, directly observable 
only for FC. Due to the disorder and the extremely weak 
anisotropy [T], the spin structure is not able to lock to 
the underlying crystalline structure for ZFC so that only 
rings of intensity can be observed as shown in Fig{l]D. 
As shown in Ref. [1 the A-Phase can be identified as a 
hexagonal lattice of skyrmion tubes. The trend to align 
two intensity maxima along (110) for field in the (111) 
and (110) directions can be easily explained in the model 
described in Ref. [1 by adding, e.g., an anisotropy term 
"^ki^x ky -\- k^)rnj^m^j^ with a positive coefficient. 

The observation of two domain populations for (100) 
(Fig.[l]G), which change as a function of field and temper- 
ature, unambiguously distinguishes the skyrmion lattice 
from a single-Q multi-domain state. However, by sym- 
metry, the contribution of the 6th order anisotropy term 
discussed above (and similar terms with a 7r/2 rotation 
symmetry) vanishes to leading order for a field parallel 
(100). Indeed, the pattern shown in FigjlJC differs qual- 
itatively from Figs{l]E and[l]F. As the Co atoms reduce 
the symmetry locally one may speculate that pinning 
terms of random sign arising from the disorder explain 
the two six- fold domain populations shown in Fig|l]G. 

Our data further allude to a possible analogy of 
skyrmion lattices with superconducting ffux line lattices 
[3j. For instance, the superconducting ffux lattice in Nb 
shows two domain populations for (100) driven by the 
frustration between the sixfold symmetry of the ffux lat- 
tice and the four-fold symmetry of crystal lattice (see 
e.g. [31 ). Likewise, the A-phase for (100) shows two do- 
main populations, thus indicating a ffat potential land- 
scape, which may also explain the peculiar regime of the 
A-phase for (100) (Fig.[2]C). Further, the metastable fea- 
tures of the A-phase are similar to superconducting ffux 
lattices in the presence of disorder. For instance, in the 
presence of disorder the superconducting ffux lattice fre- 
quently leads to a ring of intensity in SANS [32 , reminis- 
cent of the A-phase after ZFC. Similarly, for FC the su- 
perconducting ffux lattices survives as a metastable state 
in the Meissner phase. This is akin to the ffeld range in 
which the A-phase during FC is observed down to the 
lowest temperatures (we attribute the temperature de- 
pendence well below Tc to that of the ordered moment). 

An outstanding puzzle is the intensity distribution of 



the ZFC state, which is a priori not connected with 
the skyrmion lattice in the A-phase. However, two sce- 
narios have been suggested with a similar context: a 
metastable body centered multi-Q state [27J, or weakly 
stratiffed skyrmions lines [28^. Thus the ZFC state and 
its metastable nature raise two exciting issues for future 
studies. First, whether disorder enhances the stability of 
skyrmion lattices in an applied ffeld and, second, whether 
skyrmions may even form spontaneously in zero ffeld. 

In conclusion, our small angle neutron scattering study 
of Fei_a^Coa^Si identiffes a second material after MnSi in 
which a skyrmion lattice forms. As its most important 
facet, the skyrmion lattice appears in a strongly doped 
semiconductor, where disorder leads to hysteretic behav- 
ior and the observation of a skyrmion lattice in a wide 
temperature range. 
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